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ABSTRACT 

Aims. We study the supernova remnant 0506-68 in order to obtain detailed information about, among other things, the ionisation state 
and age of the ionised plasma. 

Methods. Using the Reflection Grating Spectrometer (RGS) onboard the XMM-Newton satellite we are able to take detailed spectra 
of the remnant. In addition, we use the MOS data to obtain spectral information at higher energies. 

Results. The spectrum shows signs of recombination and we derive the conditions for which the remnant and SNR in general are 
able to cool rapidly enough to become over-ionised. The elemental abundances found are mostly in agreement with the mean LMC 
abundances. Our models and calculations favour the lower age estimate mentioned in the literature of ~ 4000 year. 

Key words. ISM: supernova remnants - supernovae: general 



1. Introduction 

Supernova remnants (SNRs) hold important information about 
the nucleosynthesis and energy of the supernova explosion, its 
circumstellar matter (CSM) evolution and the physics of the 
ionised plasmas. The Large Magellanic Cloud (LMC) is partic- 
ularly well suited for the study of SNRs, since the absorption 
column to the LMC is low and it is relatively close (50kpc). 
Moreover, studying SNRs in the LMC has the advantage that the 
distance is known quite precisely, which is convenient for lumi- 
nosity and length scale calculations. 

SNR 0506-68 (also known as N23) is a small (R ~ 10 pc, see 
Fig. [TJ remnant located fairly centrally in the LMC. The X-ray 
emission of the remnant is characterised by a filamentary struc- 
ture with some bright spots. There is a gradient in brightness 
runnin g from southeast to no rthwest. The remnant has been stud- 
ied bv | Hughes et aD d2006|) with the Chandra telescope. Both 
Hug hes et al.l d2006h and lHavato et al.l (120061) report the pres- 
ence of a compact object and conclude that the SNR is a re- 
sult of a core-collapse supernova explosion (SNe). They find that 
the X-ray emission comes largely from the swept-up interstellar 
medium and estimate the age of the remnant to be ~4600 yr. 
Recently. ISomeva et alJ l2010) used the XIS instrument onboard 
the Suzaku telescope to obtain a more detailed spectrum of SNR 
0506-68. They note the presence of a cool (~0.2 keV) temper- 
ature component in addition to a hot (~0.6 keV) component in 
the plasma, and a high ionisation parameter n e t (~10 13 cm~ 3 s). 
From a Sedov analysis they conclude that the age of the rem- 
nant, based on the cool component, may be as high as ~8000 yr, 
but that it may have entered the radiative phase of its evolution. 
If the age is indeed as high as 8000 yr that would mean that the 
size of the remnant is quite small for its age, suggesting that the 
explosion took place in a high density region of the LMC. 

It has long been recognized that when a gas is suddenly 
heated in a shock it is wniier-ionised, and requires a density- 
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Fig. 1: A smoothed Chandra RGB image of the SNR 0506-68. The ar- 
rows denote the direction of the RGS dispersion directions. The 
RGS field of view covers the remnant completely. Red corre- 
sponds to Ovn (0.53-0.61 keV), green to FeL (0.79-0.89 keV) 
and blue to Mgxi (1.25-1.41 keV). The arrows denote the dis- 
persion axis orientation of the 2000 (upper arrow) and the 2002 
(lower arrow) observations. 



weighted time scale n e t ~ 10 cm s to approach equilibrium 
(Smith & Hughes 2010). Recently, however, two groups have re- 
ported evidence for over-ionised plasma in the SNRs W49B and 
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IC 443 (lYamaguchi et al]l2009t lOzawa et alJl2009l: iMiceli et ail 
l2010h . In this paper we aim to obtain detailed understanding 
about the ionised plasma of SNR 0506-68 using the spectral di- 
agnostic capac ity of the RGS (Reflection Grating Spectrometer) 
instrument dden Herder et alJfeOOll) . High resolution spectra of 
SNRs, obtained with grating spectrometers such as the RGS, 
hold detailed information about the ionised plasma. The spectral 
resolution of the RGS is large enough to resolve, among others, 
the OVII He-a triplet at ~ 22 A and the Fe xvn line complex 
around 15-17 A, provided the source has a small angular extent. 
The ratios of different emission lin es in these elemen ts provide 
interesting plasma diagnostics (e.g.: lPorquet et alJl2010h . 

2. Data 

We used the data obtained by the XXM-Newton satellite on 
July 6, 2000 (obs ID 0111130101) and July 10, 2002 (obs ID 
011130701). The exposure times of the observations are 18.1 
and 19.7 ks. Fo r our spectral study , we made use of the RGS 
and EPIC MOS dTurner et alJ200ll) data. Although the EPIC pn 
camera has a higher sensitivity, the MOS cameras have a higher 
spectral resolution, which is important for comparison with the 
RGS data. All MOS observations were performed in full frame 
mode with the medium filter in place. 

The 2000 RGS observation was taken with the dispersion di- 
rection oriented 27° counterclockwise from the celestial north, 
while the 2002 observation was rotated 90° with respect to this 
observation, at a dispersion direction orientation of 117° (see 
Fig. [TJ. This resulted in different line profiles which need sep- 
arate responses, as the 2002 observation has a higher effective 
resolution than the 2000 one. The RGS data were corrected for 
periods of high background flaring by creating good time inter- 
vals based on the count rate in CCD number 9 of the instru- 
ment. This CCD is closest to the optical axis of the telescope 
and therefore affected most by background flaring. The second 
order spectra have a higher resolution, but are of lower statis- 
tical quality. Since they provide no additional information they 
are not presented here. 

The RGS is a slitless spectrometer. When using this kind of 
instruments, the lines in the spectra get smeared out as a result 
of the extent of the source. Although the angular size of SNRs 
observed in the LMC is modest, this smearing is still present. 
We corrected for this using the heasoft program rgsrmfsmooth 
(A. Rasmussen). This program calculates a brightness profile, 
based on an image, in the direction of the RGS dispersion axis. It 
then uses this profile to adjust the response matrix to correct for 
the line smearing. In our analysis a Chandra image was used to 
create the profile, to obtain maximum precision in the correction. 

The spectra wer e analysed using the SRON SPEX package 
dKaastra et alJ [l996). Before the spectra were fit, we performed 
the SPEX optimal binning. This bins the data into optimal bins, 
making use of the statistics of the source as well as the instru- 
mental resolution. 

All standard RGS and MOS reduction tasks were done using 
XMM SAS version 10.0.0. 

3. Results 

3.1. Overall spectra 

The total RGS spectrum of SNR 0506-68 is shown in Fig.|2]and 
the MOS spectrum is shown in Fig. [3] As is clear from Fig. [2] the 
spectrum is dominated by emission lines of highly ionised O, Fe, 
Ne, N and C . The oxygen emission is particularly present with 



the notable O VII line triplet at ~ 22 A and the O VIII lines at 
~19 and ~16 A. The fact that both Ne and C are present sug- 
gests that there exis t both a cool and ho t component in the local 
plasma. Thus, like lSomeva et al.l (I2010). the spectra were fit us- 
ing two non-ionisation equilibrium (NEI) models. This model 
attempts to fit the data with two values of the parameters emis- 
sion measure n e Vx, electron temperature T e and ionisation 
age n e t. We used the standard SPEX absorption model to pro- 
duce the hydrogen column A^h to the remnant. Before fitting with 
the RGS, we first obtained the high T continuum slope by fitting 
the MOS data. When fitting the RGS data, the abu ndances were 
first fix ed at the LMC / ISM abundances found bv lHughes et al.l 
(Il998l) . and the ionisation parameters of the two NEI compo- 
nents were coupled. When the different temperatures of the mod- 
els had converged, the ionisation parameters as well as some of 
the abundances were released to obtain a better fit. 

For the RGS data, models with different ionisation 
timescales n e f and temperatures T give more or less comparable 
fits to the data. The models that were investigated can roughly 
be divided in four categories which are listed in Table Q] The 
model listed at the bottom row of the table (model D) consists of 
two NEI components which are in ionisation equilibrium, i.e: 
n e t > 10 12 cm" 3 s. This model is c omparable to the best fit 
model found by Somev a et al.l (|2010). A lower C-statistic was 
found for a two component NEI model, in which one of the com- 
ponents had a somewhat lower ionisation timescale (model C). 
This indicates that at least part of the remnant's plasma is out of 
ionisation equilibrium. The two best-fit models of table Q] are a 
two component NEI model, in which both components have a 
low n e ? (model A) and a cooling model (model B, explained be- 
low). Because the difference in C-statistic between these models 
is very minor, considering the amount of degrees of freedom, we 
listed the parameters of both the models in table [2] The ionisa- 
tio n parameters of mod el A and B are comparable to those found 
bv lHughes et ail (120061) . 

The best fit model of the 2002 RGS spectrum, the cooling 
model (model B), deserves special attention. This is a model 
in which one of the NEI components is inverted, i.e.: the ini- 
tial temperature of this component is higher than the final tem- 
perature. This cooling model, with an initial temperature of 
3.0 keV, reproduces especially the Ovu resonance to forbid- 
den line ratio better than a non-cooling two temperature NEI 
model. Furthermore it produces radiative continua in the higher 
energy parts of the spectrum, which become important in the 
MOS spectra. Physically, the cooling model corresponds to a 
plasma for which the cooling rate exceeds the recombination 
rate, which can cause overionisation. As mentioned in the in- 
troduction, rapidly cooling plasmas have been observed before 
in, among others, the mixed morphology remnants IC 443 and 
W49B. The fact that the overionisation model works well for 
SNR 0506-68 could mean that overionisation is not limited to 
mature SNRs of the mixed morphology class alone. 

The MOS spectra were also fit with different combinations 
of n e f and T, similar to the RGS data. Again, model A and model 
B gave the best fit to the data. In contrast to the RGS data, how- 
ever, the difference in best fit C-statistic between model A and B 
is significant for the MOS data, namely 100 <x. When comparing 
the fits of model A and B, it is not immediately clear where the 
difference of 1 00 cr in C-stat comes from. Distinct radiative re- 
combination edges, as have been observed in e.g. W49B, are not 
clearly visible in our spectrum. Nevertheless, a more detailed in- 
spection reveals statistically significant differences. Fig.|4]shows 
the model fits in the energy range 0.8 - 1.8 keV. The dashed 
red line shows the best fit model, model B, while the black line 
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Fig. 2: The total 2001 RGS 1 and RGS2 spectrum of SNR 0506-68 in the range 5-35 A. The best fit two-component NEI model (see text) is plotted 
in red. 



shows model A. The bottom part of the figure shows the residu- 
als of the data with model A, with the difference between model 
B and A plotted as a dashed red line. This shows that model B 
follows the overall shape of the data better than model A. The 
presence of the Fe XVII recombination continuum at 1.26 keV, 
for example, improves the fit in the 1.2-1.3 keV range by lower- 
ing the C-stat by ~ 9 cr. 



As said, the abundances were coupled between the different 
model components. However, some parameters were decoupled 
to check if there were significant differences between the hot 
and cool component. The only significant difference occurred 
when the Fe abundances were decoupled. In all cases this lead 
to a significant improvement of the fit (the C-stat/d.o.f. decreased 
to 4144/2686 for our best fit model). The Fe abundance of the 
low temperature component jumped to values of five times solar 
in case of the RGS data, while the high temperature Fe abun- 
dance decreased to ~ 0.1 solar. In addition, the temperature of 
the lower T component decreased to a value of ~ 0.14 keV. It 
is possible that there is some cool Fe present in the SNR, as 
this h as been found before in mature SNRs (e.g. Uchid aet alJ 
2009). However, at a temperature of 0.14 keV the iron emission 
increases considerably when the temperature is raised by even 
a relatively small amount (~ 0.1 keV). It is possible that the 
Fe emission requires a higher temperature than other elements 
or that there are small temperature gradients present, and that 
the model compensates for this by increasing the Fe abundance. 
Since we considered a five times solar abundance of Fe non- 
physical in SNR 0506-68, we kept the abundances between the 
two model components fixed. The RGS is less sensitive at higher 
energies, so the MOS data were used to constrain the abundances 
ofMg, Si and S. 



Table 1: RGS data C-stat / d.o.f. of the tried models for the different 
observations. 



Model Mi (high kT) n e t 2 (low kT) 2000 



2002 



A 
B 
C 
D 



low 
low 
low 
high 



low 
low 
high 
high 



4801 /3160 4228/2684 

4806/3161 4218/2685 

4820/3160 4234/2684 

4832/3159 4282/2683 



Notes. All models consist of two NEI components which had differ- 
ent kT; one high (0.6-0.85 keV), and one low (0.1-0.25 keV). The dif- 
ferent components also have different n c t, high ~ 10 14 cnT 3 s and low 
~ 10 10 ~ H cnr 3 s. Model B is a cooling model (see text). The parameters 
of models A and B are listed in Table 2. 



3.2. Detailed line spectroscopy 

As shown above, model B, i.e. a cooling plasma, gives a good fit 
to our data. If the plasma is indeed rapidly cooling, there could 
be some other spectral indications. Below we will investigate 
several known plasma diagnostics to obtain detailed information 
about the spectrum. 

An important diagnostic is the ratio of the different lines of 
the O vii line triplet. In the 2002 spectrum, the forbidden line at 
22.08 A is underestimated by the non-cooling model, while the 
resonance line at 21.6 A is overestimated. As mentioned above, 
this line triplet was best reproduced by manually tweaking an 
NEI model to mimic a recombining plasma. Recombination 
preferentially populates the triplet levels that feed the forbid- 
den line, while the resonance line is populated mainly by col- 
lisional excitation. As such, an enhanced forbidden line suggests 
the presence of enhanced recombination in the plasma. This will 
be discussed in more detail in paragraph s. 31 
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Model A Model B 



Parameter 






RGS 


MOS 


RGS 




MOS 


1 * H 


(10 21 cm" 


I 


1.14 ± 0.01 


1.14 ± 0.01 


1.14 


±0.01 


1.14 ± 0.01 


n^uu V\ 

,t £' t Ii ' 1 


(10 58 cm" 


-3) 


20.5 ± 0.5 


14.9 ±0.2 


10.5 


±0.3 


1 1.4 ± 0.1 


n e n h V"> 


(10 58 cm" 




99.6 ± 3.7 


103.5 ± 1.6 


190: 


t7.5 


138.5 ± 1.6 


Pre shock kT 


(keV) 








3.0 




3.0 


kT\ 


(keV) 




0.85 (fixed) 


0.85 ±0.01 


0.85 (fixed) 


0.85 ± 0.01 


n £ t[ 


(10 10 cm 


" 3 s) 


4.67 ± 0.32 


7 20°- 06 

'•^"-0.17 


? 3? 


+ 0.13 


6.7 ±0.2 


kT 2 


(keV) 




0.19 + 0.00 


0.22 ± 0.00 


0.15 


±0.00 


0.18 ±0.00 




(10 10 cm 


" 3 s) 


41.2 + 6.6 


52 l 610 

_ 39 


54.5 


±3.7 


99.6 ± 2.7 


L x 


(10 34 erg 


s- 1 ) 


12.5 


9.4 


14 




8.8 


C-stat / d.o.f. fl 








1170/632 






1068/633 


Element 






Abundance (wrt solar) 








V 






0.55 ± 0.08 




0.18 


±0.03 




N 






0.06 ± 0.02 




0.07 


±0.02 











0.20 ±0.01 


0.19 ±0.01 


0.26 


±0.01 


0.22 ± 0.01 


Nc 






0.22 ± 0.02 


0.22 ± 0.02 


0.33 


±0.03 


0.24 ± 0.02 


Mg 








0.31 ±0.02 






0.27 ± 0.02 


Si 








0.25 ± 0.03 






0.19 ±0.02 


S 








0.33 ± 0.09 






0.31 ±0.07 


Fc 






0.26 ±0.01 


0.29 ±0.01 


0.38 


±0.02 


0.43 ± 0.02 



Notes. The Anders & Grevesse ( 1989) Solar abundances were used.The abundances are comparable between the instruments. 
(fl) The C-stat/d.o.f. of the RGS data can be found in Table 1. 




Energy (keV) 



LU 

"CD 



0.1 



S 0.4 



^3= *± 






Model A 

Model B - 








j T 



0.4 I 1 1 _ 1 1 — t 



1.2 1.4 

Energy (keV) 



1.6 



Fig. 3: Plot of the MOS 1 and MOS 2 data of 2002. The model seems 
to fit the data well. 



All tested models have trouble fitting the FeXVll 15-17 A 
line ratio (see Fig. 01. These lines are formed by the transi- 
tions from the 3d and 3s levels of the Ne-like ion to the ground 
state. We used the same labeling of the lines that was used 
bv iGillaspv et all (1201 ll) . The fact that the 17 A line blend is 
stronger than the 15 A line can b e another sign of enhanced re- 
combination dLiedahl et al.lll990h . To obtain the Fe xvn line ra- 
tios, the best fit overall model was used, while the contribution to 
the emission by Fe xvn lines was excluded from this model; this 
ensures that contributions from, e.g., higher lines of the O VIII 
Ly series are taken into account. The Fe xvn line complex in the 
range of 14-18 A was then fit with five gaussians, with a fixed 



Fig. 4: A plot of the mos 2 data with model A and model B in the 
range 0.8 - 1.8 keV. The bottom part of the figure shows the rel- 
ative error between model A and the data, while the dashed red 
line shows model B - model A. It is clear that model B follows 
the overall shape of the data better than model A. For example, 
model A shows a residual at ~ 1.25 keV, whereas model B im- 
proves the fit of this region due to the Fexvn recombination 
edge at 1.26 keV. 



Nu of 1.14 x 10 21 cirT 2 . The observed line strengths are listed 
in table [3] We can compare our Fe xvn line ratios with recent 
laborator y measurements to obtain more information about the 
plasma. Gillas py et al.l(l201 ll) measure the 3s/C and C/D ratios at 
different electron beam temperatures, where 3s = F+G+H. Our 
observed 3s/C ratio of 2.45 + 0.2 and C/D ratio of 1.46 + 0.13 
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0.3 



It 



14 14.5 15 15.5 16 16.5 17 17.5 

Wavelength (A) 

Fig. 5: The 2002 RGS1 and RGS2 spectrum of SNR 0506-68 in the 
range 14-18 A. The best fit two-component NEI model (see 
text) is plotted in solid black. The model has trouble fitting the 
15 and 17 A Fexvn as well as the O vm Ly pfLy a line ratio. 

Table 3: Line fluxes of a number of emission lines calculated with and 
without hydrogen absorption. 



Notes. (a) In units of 10 44 ph s~ 




0.0008 0.001 0.0012 

Hydrogen absorption column (10 24 cm" 2 ) 



Fig. 6: Contour plot of the carbon abundance and the Hydrogen absorp- 
tion column made using model B. The contours represent the 1 
and 2 a confidence regions. 



line 


A ( A) 


Intrinsic Flux" 


Flux without absorption" 




Fe xvn C 


15.01 


1.11 ±0.09 


0.76 ± 0.06 




Fe xvn D 


15.24 


0.76 ± 0.09 


0.50 ± 0.06 




Fe XVII F 


16.78 


0.52 ± 0.08 


0.32 ± 0.05 




FexvuG// 


17.05 


2.21 ±0.16 


1.22 ± 0.09 




O vm Ly/? 


16.01 


1.00 ±0.09 


0.65 ± 0.06 




O vm Ly a 


18.97 


7.63 ± 0.25 


3.88 ±0.13 


Z3 


Ovilr 


21.60 


5.90 ± 0.41 


2.30 ±0.16 


O 

o 


O vil i 


21.80 


1.37 ±0.28 


0.55 ±0.10 


X 
13 


Ovu/ 


22.10 


4.69 ± 0.38 


1.69 ±0.14 


Li- 


CviLyyS 


28.47 


0.48 ±0.19 


0.19 ±0.06 




C VI Ly a 


33.74 


5.83 ± 0.80 


0.71 ±0.10 





0.15 - 




0.05 - 



21.4 21.6 21.8 22 
Wavelength (A) 



21.4 21.6 21.8 22 
Wavelength (A) 



both correspond to an electron beam temperature of ~ 0.85 keV, 
which is near the collisional excitation threshold. The observed 
line ratios show no indication of a high 3s/3d ratio indicative of 
recombination. 

At 16 A, the O vm Ly f3 line is underestimated by the model. 
By fitting the 15-20 A region with an absorbed continuum and 
gaussians, a Ly /6/Ly a ratio of 0.13+0.01 is obtained. This value 
is in agree ment with CIE val ues found for this ratio at T - 4 - 
5 x 10 6 K (ISmith et aUl200lh . 

C vi at 33.7 A also has a strong presence and is detected at 
a 9cr level. Because it is such a strong line, and it is affected 
heavily by absorption, it can be used to constrain the absorption 
column to the remnant. This was done by making a contour plot 
(Pig. O of the Nh and the C abundance after a good fit was ob- 
tained. This resulted in an Nh of 1.14 x 10 21 cirT 2 , which was 
used in all our models. Note that in our model the C vi originates 
from the coolest component only (kT = 0. 15-0.20 keV). At these 
temperatures C is mostly ionized, with the C vi fraction being as 
low as 25%. In principle additional Cvi emission could come 
from an even lower temperature component. However, since the 
cooling timescale for such a component is very short a major 
contribution does not seem very likely. 



Fig. 7: The 2000 and 2002 observation of the OVII He-a triplet. It is 
clear that the 2002 observation has a higher effective spectral 
resolution. 



3.3. G-ratio 

Fig. [7] shows the 2000 and 2002 observations of the Ovu He-a 
triplet. This triplet consists of a resonance, forbidden and inter- 
combination line. F (A — 22.098) li2i 3 5 1 — > Ls 2 l So is the for- 
bidden transition, I (A = 21.804, 21.801) is the sum of the two 
intercombination transitions \s2p i P\2 — > Is 2 l So, and R (A = 
21.602) is the resonance transition ls2p l Pi — > Is 2 1 Sq. Both 
the 2000 and the 2002 observation are well-fit by three gaus- 
sians with a fixed Nh- There are some differences between the 
two observations, the most notable one being that the F/R ra- 
tio is larger in the 2002 observation. In addition, the lines in the 
2000 observation are broader, which is an effect of the orienta- 
tion of the RGS dispersion axis. The total flux in the line triplet 
is approximately equal between the observations. 

An interesting quantity which can be de rived from these line s 
is the so-called G-ratio; G = (F+I)/R (e.g.: lPorquet et al.l2 010). 
This quantity equals 0.87 + 0.09 for 2000 and 1.19 + 0.09 for 
2002, giving a combined ratio of 0.99 + 0.06. As the 2002 ob- 
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servation has a higher effective resolution, the G-ratio of that ob- 
servation may be more reliable. Fig. [8] shows a plot of G-ratios 
for different n e f. The values of the G-ratio for different tempera- 
tures reach a constant value as the plasma approaches collisional 
ionisation equilibrium (CIE). If we take the mean value of the 
above G-ratios, the calculated G-ratio at high n e t and kT = 0.2 
keV lies just within the error bars and is thus as expected. If 
the 2002 observation is indeed more reliable, the deduced G- 
ratio of the plasma lies above this CIE value and the plasma is 
over-ionised and could be recombining. In principle the excess 
photons present in the forbidden line should show up as a recom- 
bination edge at 16.78 A. A recombination edge in addition to 
the already present continuum was not found, however. 



4. Discussion 

We made a detailed spectral analysis of the SNR 0506-68 using 
mainly the RGS instrument aboard the XMM-Newton telescope. 
The best fit to the overall mos and RGS spectra is model B: a two 
component NEI model, of which one component is inverted. We 
investigated the hypothesis that the plasma is cooling leading 
to an overionisation, using some known plasma diagnostics. Of 
those diagnostics, only the 2002 O vn triplet line ratio confirms 
the hypothesis. There are however, several other physical mech- 
anisms which can cause the observed line ratio. 

Resonance scattering can cause photons of resonance lines to 
be scattered in the direction of least optical depth, reducing the 
line flux if the optical depth along our line of sight is high. This 
process was suggested to be respo nsible for the observed O vn 
F/R ratio in the SNR DEM L71 (I van der Hevden et all 120031) . 
The scattered photons are not lost, however, so the remnant must 
have a specific geometry, i.e. it cannot be spherically symmetric, 
or resonance scattering will have no effect whatsoever. The op- 
tical depth in the O vn line equals about 6 in our line of sight, 
without taking microturbulent velocity into account. The opti- 
cal depth equals 1, however, at a turbulent velocity of 80 km/s 
and decreases even more at higher values. A high optical depth 
could result in a signific ant reduction of the resonance line flux 
dKaastra & Mewd fl995). which means that resonance scattering 
could be significant in this remnant. 

Charge exchange occurs when a highly ionised gas enters a 
neutral gas region, such as a molecular cloud. H-like ions collide 
with the neutral ions to form excited He-like ions which decay 
by a radiative cascade. This process enhances the forbidden line 
due to the higher statistical weight. There have been no reports 
of neutral gas regions close to the remnant, so the probability of 
the charge exchange scenario occurring remains uncertain. 

4.1. Cooling 

As mentioned above, recombination can influence the observed 
line ratios. For recombination to dominate, the plasma must be 
over-ionised and thus must have cooled faster than the recombi- 
nation process happens. We now investigate if SNR 0506-68 can 
cool fast enough for this to occur. 

The expanding plasma cools in two ways: through radiation 
and due to adiabatic expansion. A timescale for radiative cooling 
can be obtained by dividing the internal energy of a gas by the 
luminosity: 
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Fig. 8: G-ratio of O vn at different temperatures, as a function of ioni- 
sation parameter. These curves where calculated using SPEX. 



where A is the cooling rate in units of erg cm 3 s _1 , k is 
Boltzmann's constant, and «h is the hydrogen/proton number 
density. 

For an adiabatic gas, PV r = constant. By the same token, 
TV 7-1 = constant. The following formula for the adiabatic cool- 
ing timescale can thus be derived for a spherical volume: 



Tad 



V r/ad \r - 1/ 3/? 



(2) 



A total cooling timescale as a result of radiative and adiabatic 
cooling is now given by r"^ = r~l + r~J . When a plasma cools, 
the radiative cooling rate of the remnant is proportional to the 
density squared. The recombination rate also increases propor- 
tional to the density squared. If adiabatic cooling is not important 
compared to radiative cooling, the remnant will stay in ionisation 
equilibrium. Radiative cooling is faster than recombination be- 
low 10 6 K or if the elemental abundances are strongly enhanced. 
For an over-ionised plasma, it is required that r rec > r coo i. 
The recombination timescale for O vin — > O vn in the range 
10 6 < T < 10 7 is approximately equal to (1.3 x 10 _12 n e ) _1 - 
Using [T|[2] with y = 5/3, this leads to: 



1 



2R nA 
Y + 3kf < 1.3 x 10- 12 n e 



(3) 



Trad 



«hA 



(1) 



Note that the remnant needs to have reached ionisation equi- 
librium before adiabatic cooling causes an over-ionised plasma. 
For/? =* 10 pc, R ^ 400 km s~\ A = 3xl0~ 23 erg c m 3 s" 1 , corre- 
spond ing to a 0.2 keV plasma at LMC abundances dSchure et al.l 
2009) and a density n — 10 cm -3 , the recombination timescale 
(using the above formula) is about five times lower than the cool- 
ing timescale, and depends strongly on the density. 

To expand the above equation to a larger temperature range, 
we can rewrite R/R = fi~ x t, where /3 is the expansion parameter 
and t is the age of the remnant. Values of (3 vary between about 
0.4 for a remnant in the Sedov expansion phase and 0.25 in the 
snowplough phase. In addition, we can use the temperature de- 
pendent O viil — > O VII (radiative + diel ectronic) recombination 
rates from Shull & van S teenberg ( 1982) as w ell as th e complete 
O radiative cooling curves from lSchure et al.l (12009b . 

Fig. [9] shows this relation between the cooling and the re- 
combination time using t - 10 11 s and ji = 0.4. The behaviour 



6 



S. Broersen et al.: The high-resolution X-ray spectrum of SNR 0506-68 using XMM-Newton. 



1000 



100 



10 



n=1 
n=10 
n=100 



0.1 



5.5 



6.5 
Log T 



fore estimate the volume of the sphere to be 9 x 10 58 cm 3 . The 
fact that the radiation is distributed anisotropically over the rem- 
nant as well as the fact that the emission does not come from a 
purely spherical region, can be accounted for by using a filling 
factor f. A filling factor within 10-25% seems reasonable. 

Using the normalisation factor n e rtuVx of the best fit model 
to the RGS spectra, with « e = «h /1.2 and Vx the volume, we 
obtain a density «h = H(/o.i) 0,5 cirr 3 and «h = 5(/o.i)~ ' 5 
cm -3 for respectively the cool and the hot component of the 
non-cooling NEI model, where /o.i represents a filling factor 
of 10%. The density of the cool component of the cooling 
model is slightly higher at 14(/o.i) ■ The density of the low 
T component is som ewhat higher than previously found values. 
Hug hes et al.l (120061) find de nsities of 10-23 cirr 3 using X-ray 
data. iDickel & Milnd (1 19981) used radio data to find a maximum 
density of 10 cirr 3 and lWilliams et al.l (120061) found a density of 



Fig. 9: A comparison of the oxygen recombination timescale with the 
cooling timescale at different densities. Low density favours the 
appearance of an overionised plasma, because the recombination 
timescale lowers proportional to the density, as does radiative 
cooling, but radiative cooling only becomes important at log T < 
10 6 K. These curves where calculated using LMC abundances. It 
should be noted that a plasma at LMC abundances is less likely 
to reach an over-ionised state at T < 10 6 K, because of the lower 
radiative cooling rate. 



at high temperature is dominated by the adiabatic cooling, while 
radiative cooling becomes important at T < 10 6 K. An increased 
density increases the recombination rate, which decreases the 
likelihood of plasma in SNRs to become over-ionised at tem- 
peratures at which radiative cooling is unimportant. It should be 
noted that the above derivation is a first order estimate. We can 
however still use it to check if the conditions in SNR 0506-68 
are likely to cause an over-ionised plasma. At a T ~ 10 6 K and 
n ~ 10 cirr 3 the cooling timescale is about two times lower 
than the recombination timescale. This is a discrepancy with the 
above calculation, so we can conclude that either our remnant 
is not in the Sedov stage of evolution (i.e. (3 < 0.4), or that the 
shock velocity or density estimations are not right. 

We can check if the overionisation observed in other SNRs 
can be explained using the above relation. In W49B, recombin- 
ing Fe xxv - xxvi was found. Using their respective recombi- 
nation and radiative cooling rates at T = 1.5 keV and n — 10 
cm -3 , the cooling timescale is indeed lower than the recombina- 
tion timescale. 

Note that the above derivation for the cooling rate is some- 
what conservative, as we did not take into account dust cool- 
ing. This is in parti cular true for SNR 0506-68, for which 
Williams et al. (2006) found an infrared luminosity of 8.7 x 10 36 
erg s . This is about 60 times higher than the observed X-ray 
luminosity. In general dust cooling increases the cooling rate of 
the plasma, which makes it easier to reach an over-ionised state. 



4.2. Explosion parameters 

The data obtained from fitting allow us to constrain several in- 
teresting parameters. For this we need the volume of the plasma, 
which we obtained using the Chandra data. The radius of the 
sphere in which the plasma is contained is about 38", which cor- 
responds to a radius of 9.1 pc at a distance of 50kpc. We there- 



7-5 8 5 8 cirr 3 by modelling dust emission in the remnant. These are 
mean densities over the whole remnant and it is probable that at 
some locations the densities are higher. 

With these densities and the ionisation parameters the shock 
ages of the cool and hot components are ~ 1200 and ~ 115 
yr. Different parts of the remnant were thus shocked at different 
times, which suggests that there are some density fluctuations 
in the surrounding IS M, which is also clea r from Fig. 1. It has 
already been hinted by Hugh es~et al.l d2006) that the open cluster 
HS 114 (Hodg e & Sextonl 19661) . which lies on the brighter, high 
density side of the remnant, may be the cause of the observed 
brightness gradient. 

The total swept up mass is given by M swept ~ num p fV, 

where m p is the proton mass. M swept ~ 88 -J/o.i M Q for the cool 



component and M s 



wept 



37 y/o.i M for the hot component. 



Due to the brightness gradient and the fact that a low and 
high temperature NEI model are necessary to fit the spectrum of 
the remnant, an age estimation is difficult. Using a Sedov model 

(fsedov = 4.3 x 10 2 (j^) ( jf^) ' ), the age of the cool compo- 
nent can be estimated at ~ 9000 year, while the age of the hot 
component equals ~ 4000 year. A Sedov model assumes that the 
remnant is expanding in a homogeneous medium which, judg- 
ing from the anisotropic emission, is not quite valid in the case 
of SNR 0506-68. Since the hot, less dense component of the 
plasma is likely to have expanded relatively undisturbed, the age 
estimation of this component is probably more accurate. 

5. Conclusions 

- The SNR 0506-68 is best fit by a two component NEI model, 
of which one of the components is cooling. This suggests 
that the SNR may be cooling, but the emission line diagnos- 
tics marginally agree with this scenario. 

- The most natural explanation for the enhanced G-ratio of 
the O vii He-a is that part of the plasma is recombining, al- 
though resonance scattering cannot be ruled out. 

- When an expanding plasma reaches near-ionisation equilib- 
rium, adiabatic cooling can cause it to become over-ionised. 

- The abundances of SNR 0506-68 are very similar to the 
mean LMC abundances. This, coupled with the fact that 
Mswept is fairly high, confirms that the emission is dominated 
by the emission from the shock heated ISM. 

- Overionisation may be more widespread in mature SNRs 
than usually thought. 

- The age of the remnant is uncertain and model dependent. 
However our models and calculations favour the lower age 
estimate mentioned in the literature, namely ~ 4000 year. 
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